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<N 1 Introduction 

A perfect matching of a graph is a subset of edges which covers every vertex exactly once, 
O ' that is, for every vertex there is exactly one edge in the set with that vertex as endpoint. 

The dimer model is the study of the set of perfect matchings of a (possibly infinite) graph. 
The most well-known example is when the graph is Z 2 , for which perfect matchings are 
equivalent (via a simple duality) to domino tilings, that is, tilings of the plane with 2x1 
and 1x2 rectangles. 

i ' In the first three sections we study domino tilings of the plane and of finite polygonal 

^ ■ regions, or equivalently, perfect matchings on Z 2 and subgraphs of Z 2 . 

In the last two sections we study the FK-percolation model and the dimer model on a 
more general family of planar graphs. 

y—i ! 

m ■ 

2 The number of domino tilings of a chessboard 

^ ■ 

A famous result of Kasteleyn [S] and Temperley and Fisher ^H] counts the number of domino 
tilings of a chessboard (or any other rectangular region). In this section we explain Kaste- 
j> ! leyn's proof. 

•i-H . 

X: 

2.1 Combinatorics 

Let R be the region bounded by a simple closed polygonal curve in Z 2 . A domino tiling of R 
corresponds to a perfect matching of G, the dual graph of R: G has a vertex for each lattice 
square in R, with two vertices adjacent if and only if the corresponding lattice squares share 
an edge. 

Theorem 1 [Kasteleyn, 1961] The number of domino tilings of R is v/j det K\, where K is 
the weighted adjacency matrix of the graph G, with horizontal edges weighted 1 and vertical 
edges weighted i = a/— T- 
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Figure 1: Random tiling of a square 



For example for the 2x3 region in Figure 12 the matrix K is 

/ i 1 \ 

1 i 1 

1 i 

1 1 ' 
1 i 1 

\ 1 i / 

whose determinant has absolute value 9 = 3 2 . 




Figure 2: 2 x 3 rectangle and dual graph. 

Proof: Since the graph G is bipartite (one can color vertices black and white so that black 
vertices are only adjacent to white vertices and vice versa) K can be written K = 



f ° A ) 

\A l J 
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so we must evaluate 

det A = ^ sgnak la(1) . . . k na(n) 

Each term in this sum is zero unless black vertex o~(i) is adjacent to white vertex % for each 
i. Therefore we have exactly one term for each matching of G. It suffices to show that two 
nonzero terms have the same sign. 

Suppose we have two matchings; draw them one on top of the other. Their superposition 
is a union of (doubled) edges and disjoint cycles. One can transform the first matching into 
the second by "rotating" each cycle in turn. In particular it suffices to show that if two 
matchings differ only along a single cycle, the corresponding terms in the determinant have 
the same sign. 

This is supplied by the following lemma, which is easily proved by induction. □ 



Lemma 1 Let C = {vo, ■ ■ ■ , i>2k-i, V2k = ^0} be a simple closed curve in Z 2 . Let mi be 
the product of the weights on edges VqV\, i>2t>3, • • • and 1712 the product of the weights on the 
remaining edges i>ii>2, . • • (recall that vertical edges are weighted i and horizontal edges 1). 
Then m 1 = (— l) n+k+1 m 2 , where n is the number of vertices strictly enclosed by c and k is 
1/2 of the length of c. 



2.2 Rectangles 

Suppose the graph G is an m by n rectangle, with vertices {1,2,..., m} x {1,2,..., n}. To 
compute the determinant of its adjacency matrix, we compute its eigenvectors and their 
eigenvalues. 

For j,k G Z define 

. irjx . nky 
f{x,y) = sin — -sin 



m + 1 n + 1 

It is easy to check that / is an eigenvector of K, with eigenvalue 

2 cos h Li cos ■ 



m + 1 ra+l 

We didn't pull this formula out of a hat: since G is a "product" of two line graphs, the 
eigenvectors are the product of the eigenvectors of the line graphs individually. 

As j, k vary over integers in [1, m] x [1, n], the / form an orthogonal basis of eigenvectors. 
Therefore the determinant of K is 

m n ■ 1^ 

det K = J I j J 2 cos — + 2/ cos 



m + 1 " Ti + 1 

j=l k=l 



Evaluating yf\ det K\ for m = n = 8 gives 12988816 tilings of a chessboard. 

One may show that for m,n both large the number of tilings is exp(Gmn/ii + 0{m + n)), 
± + ± _ JL 

32 ^ 52 72 



where G — 1 — 4z + ik — ^ + • • ■ is Catalan's constant. 
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2.3 Tori 



It will be useful shortly to compute the number of tilings for a graph on a torus as well (e.g. 
a rectangle with its opposite sides identified). Kasteleyn showed that this could be accom- 
plished by a linear combination of four determinants. In essence these four determinants 
correspond to the four inequivalent "discrete spin structures" on the graph G embedded on 
the torus. Rather than go into details, we just note that the above proof fails on a torus since 
two matchings may differ on a loop which goes around the torus (is not null-homologous) . 
The change in sign in the determinant from one matching to the next is then a function of 
the homology class in H±(T, Z/2Z) of the loop. The result is that (for m,n both even) the 
number of tilings of an m x n torus is 

Zm,n — -{ — PoO + Pal + -PlO + -Pll)i 

where 

p*r= n n ( z+ ~ +iw+ ^)- 

Note that actually P o — in this case. For details see jH] or 
One can show that 

1 G 
lim logZ mj „ = — . 

m,rwoo mn 7T 

2.4 Inverse Kasteleyn matrix 

We are primarily interested in the uniform measure on perfect matchings of a graph G. 
A very useful consequence of Kasteleyn's counting theorem is 

Theorem 2 ( [10J ) Let T = {(u>i, &i), . . . , (w^, h}~)} be a subset of dinners, (with the jth 
dimer covering white vertex Wj and black vertex bj). The probability that all dinners in T 
appear in a uniformly chosen matching is 

Pr(T) = |det(iir- 1 K6 i ))a< ij -< fc |. 

The advantage of this result is that the determinant for the joint probability of k dimers is 
of size k independently of the size of the matrix K. In particular to compute the probability 
that a single edge appears, one just needs a single element of K~ x . 

If we can compute asymptotics of K~ l for large rectangles we can then understand easily 
the limiting measures on tilings of the plane. We'll do this later. 

A more general class of measures on perfect matchings uses graphs with weighted edges. 
If edges are weighted u(e) > 0, the weight of a perfect matching is the product of its edge 
weights, and the probability measure we are interested in is that which gives a matching a 
probability proportional to its weight (for finite graphs). There is a simple generalization of 
Theorems [T] and El for weighted graphs. 

We should also point out that versions of Theorems ^ and El hold for arbitrary planar 
graphs, see and [TU] . 
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3 The arctic circle phenomenon 



The order-n Aztec diamond is the region shown. It was defined and studied in [HI ISJ- A 
random tiling is also shown. Why does it not look homogeneous, as does the tiling of a 
square? 
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Figure 3: Random tiling of Aztec diamond 

Near the corners, you only see tiles of one "type". The way to understand this phe- 
nomenon is via the height function. 

3.1 heights 

A domino tiling can be thought of as an interface in 2 + 1 dimensions. Specifically, to a 
tiling we associate an integer-valued function h on the vertices of the region tiled, as follows. 
Rather than define the height, we define the height difference between two adjacent lattice 
points. If edge vw is not crossed by a domino, the height difference h(w)—h(v) is 1 if the face 
to the left of vw is black, and —1 otherwise. If the edge is crossed by a domino, h(w) — h(v) 
is —3 or +3 according to whether the face to the left is black or white. This defines a "one- 
form", that is a function on the oriented edges of the graph G satisfying /(— e) = —/(e). 
Moreover the one-form is closed: the sum around any oriented cycle is zero (since the sum 
around any face is zero: a domino crosses exactly one edge of each face). Therefore there is 
a function h, well defined up to an additive constant, with these edge-differences. 

The height function of a tiling gives a map from Z 2 to Z. One of the principal motivations 
for studying dominos is to study the model of random maps of this sort. 
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Figure 4: Heights in a domino tiling. 

Note that a tiling determines a height function, and vice versa. Also, the height function 
on the boundary of the region to be tiled is independent of the tiling. 

For the In x In square, the boundary function is asymptotically flat, in fact assuming 
the value at a corner is 0, it alternates on edge edge between and 1 or and —1: see Figure 
13] However for the Aztec diamond of order n, it is linear on each edge: 

For a non-planar region such as a torus, the height might not be well-defined: it can 
be well-defined locally, but on a path winding around the torus it may have a non-trivial 
period. Indeed, for a graph on the m xn torus the height will have a well-defined horizontal 
period h x and vertical period h y , both integers, such that the height change on any closed 
path winding once horizontally around the torus is h x , and similary h y for vertically-winding 
paths. 

The height of a random tiling on a planar region is a random function with those given 
boundary values. It is easy to guess that the height at an interior point of a random tiling 
of a square of side n is o(n). However for the Aztec diamond it is not at all clear what the 
height of a random tiling will be. In fact it will typically lie near (within o(n)) of its mean 
value. Moreover its mean value can be determined by solving a certain variational principle 
(minimizing a certain energy). That is the content of this section. 

3.2 CKP theorem 

The theorem of Cohn-Kenyon-Propp states this in general: given a Jordan domain U with 
function u : dU — ► K, and a sequence of domains in eZ 2 converging to U, with boundary 
height functions h t converging to u, that is eh e —>■ u, then with probability tending to 1 the 
height of a random tiling will lie near (within o(n)) its mean value, and the mean value ho 
is the unique function / maximizing a certain "entropy" 

Ent(/) = / ent(d x f,d y f)dxdy. 
Ju 
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Figure 5: Heights on boundary of an Aztec diamond. 



ent(s,t) = - {L(-Kp a ) + L(-np b ) + L(np c ) + L(-Kp d )) 

7T 



where L(x) = — J* log2sint<it is Lobachevsky's function, and p a ,Pb,Pc,Pd are certain prob- 
abilities which are determined by the following equations. 



2(Pd-Pc) 

Pa+Pb+Pc+ Pd 

sin 7ip a sin npb 



s 
t 
1 

sin np c sin 7rpd- 



The proof of the first part of this theorem is somewhat standard. The harder part is the 
computation of the function ent. 

The idea of the proof of the first part is the following. Each height function is a Lipschitz 
function with Lipschitz constant 3, since \f(xx,yi) — f(x\ + l,yi)\ < 3 and similarly for the 
y-direction. In fact we have a stronger condition that when x\ — X2,yi — y<i are even, 

- 7(^2,2/2)1 < 2max{|x 1 - x 2 \, \yi - y 2 \} (1) 

Let L be the class of functions satisfying When we scale the x— and ^/-coordinates by 
e, then ef is still in L. Note that L is compact (in the uniform topology). 

On this space, the functional Ent is upper semicontinuous: its value can only increase at 
a limit point. This follows from the strict concavity of the function ent(s, t). Concavity also 
proves unicity of the maximizing function. 
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Since L is compact, cover it with a finite number of metric balls Bf{8) of radius 5. 
We will show that for any / 6 L, the number of tilings for which eh lies within 5 of / is 
exp(^-Ent(/)(l + o(l))). This determines the size of each Bf{$). Therefore when e is small, 
the number of tilings whose normalized height function lies close to / max overwhelmingly 
dominates the number of other tilings, even combined. 

Thus almost all tilings lie near / max . 

To count the number of tilings whose height function lies near a given function / e L, 
triangulate U with a fine mesh, large on the scale of e but with mesh size tending to zero. 
Since a Lipschitz function is differentiable almost everywhere, on almost all triangles the 
function / is almost linear. A tiling which lies close to / will lie close to a linear function on 
most triangles. So it suffices to count the number of tilings whose height function is close to 
a linear function. 

However a short technical lemma (essentially subadditivity) proves that if R is a trian- 
gular region (or other nice region), the number of tilings with height function lying close to 
a fixed linear function i on that region is exp(A ■ ent(£) • (1 + o(l))), where A is the area of 
R and ent(f) is a function depending only on the slope of £. 

The explicit computation of ent(s, t) comes from computing the number of tilings on a 
torus with a given slope (s, t), that is, a given horizontal height change [sm] going horizontally 
around the torus, and [tn] vertically. 

To compute this we use staggered weights a, b, c, d on the edges of Z 2 , as in Figure EH 



-o 



6- 



-o- 
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a w b 

Figure 6: Weights on torus graph 

If we put these weights in the Kasteleyn matrix then its determinant is the sum over 
tilings of the weight of a tiling (product of the edge weights). That is, y] det K\ is the 
partition function Z(a, b, c, d) of tilings with these weights. For a tiling of the torus with 
these weights, the number of dominos of type a, b, c, d determines the height change around 
the torus: the horizontal height change is (Nd — N c )/n and the vertical height change is 
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(N a — Nb)/m. Therefore once we know 



Z(a,b,c,d)= ^ Na b Nb c Nc d N " 

matchings 

we have 

Z(a, 1/a, c, 1/c) = a mh y C - nh * 

matchings 

and we need only extract the appropriate coefficient to get the number of tilings of given 
slope. 

In 0] the asymtotic formula is given 

log Z(a, b, c, d) = \ / / log (a + bz + cw + dzw) — — , 

(2m) 2 Jgi Jgi z w 

from which the above formula for ent follows. 
3.3 flats 

For the Aztec diamond, the function / max is not analytic, only piecewise analytic. In the 
four regions outside of the inscribed circle / max is linear (as you can approximately see in the 
example tiling). Inside it is analytic. On the boundary of the inscribed circle (the arctic 
circle, since it separates the "frozen region" from the "temperate" zone) the function is C l 
but not C 2 . In fact it is C 1 ' 5 except at the four boundary-edge midpoints. 

These frozen regions result from the degeneration of the ellipticity of the PDE defining 
/max at the boundary of the domain of definition. 

It is very interesting to study what happens exactly at the boundary. We won't discuss 
this here though. 



4 Conformal invariance 

The scaling limit (limit when the lattice spacing tends to zero) of the height function on 
domino tilings tends to a nice conformally invariant continuous process, the "massless Gaus- 
sian free field", a sort of two-dimensional version of Brownian motion. In this section we 
discuss the ideas behind the proof. The original references are jTIJ H2] . 

4.1 A' as Dirac operator 

We had a lot of choice in definition of the Kasteleyn matrix. If you go back to the proof, you 
will see that we could have chosen any edge weights with the property that the weights on a 
lattice square a, j3, 7, 5 satisfy the condition that 07 / '(38 is real and negative. It is convenient 
to put weights as in Figure [7| 

Lemma ^ holds for these weights as well, as do Theorems and EJ 
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Figure 7: 



Let us now study K~ x for these weights. Because of bipartiteness, K~ 1 (vi,v 2 ) = unless 
one of Vi,V2 is black and the other is white. Let B denote the black vertices and W the 
white vertices. Let w G W . If / G C B , then Kf = means 



that is u d x + idy'f = 0. We say / G C B is discrete analytic if Kf = 0. 

Let B , Bi be the black vertices whose coordinates are both even (respectively, both odd). 
Similarly let W , W\ be the white vertices whose coordinates are (1,0) mod 2, respectively 
(0, 1) mod 2. 

Lemma 2 If f G C B is discrete analytic, then f is harmonic on Bq and B\ separately, that 



Moreover, these two harmonic functions are conjugate. 

The proof follows from the observation that K*K is the Laplacian when restricted to any 
of the sublattices B ,Bi, W , W\, see Figure |H1 

As a consequence K~ l can be related to the Green's function. 

4.2 K and the Green's function 

For a fixed w, K~ l (w,b) (considered as a function of b) is discrete analytic, except at w. 
This follows from KK^ 1 = I. 

In fact K~ x can be written in terms of the Green's function G. Recall that G(u,v) is 
defined to be the function satisfying AG(u,v) = 5 u (v) (taking the Laplacian wrt the second 
variable). For a fixed w, 



f(w + 1) - f(w - 1) + i(f(w + i) - f(w - i)) = 



IS, 



4/(6) = fib + 2) + fib - 2) + f{b + 2i) + fib - 2i). 



AK-\w, ■) = K*KK~ 1 (w, •) = {K*5 W ){-) = 8, 



'w+l — Ow-1 ~ 
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Figure 8: 
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Thus K 1 (w, •) is a sum of four Green's functions. When restricted to B , it is the difference 
of two Greens functions, K~ 1 (w, 6) = G(w + 1, 6) — G(w — 1,6). 

4.3 On a rectangle 

We can easily compute K~ x for a rectangle, since we already diagonalized the matrix K. 
Because K is symmetric, we have the following: 

K- 1 (v 1 ,v 2 ) = Y,fMfM/X j 



where fj are the orthonormalized eigenvectors with eigenvalue Xj. 

In the limit when the rectangle gets large, and t>i, v 2 are far from the boundaries, we can 
find 

i p2ir i-2-k i(xe+y<j>) 

K-\v u v 2 ) = — / — — - — d9dc/>, (2) 



4tt 2 Jo Jo 2z sin 9 — 2 sin < 

where (x, y) = v 2 — v\ G Z 2 . 

This formula defines, along with Theorem El a measure on tilings of the plane. This 
measure is the unique weak limit of measures on rectangles. It is also the unique entropy- 
maximizing translation-invariant measure on tilings [2]. Some values of K~ l are shown in 
Figure 

Asymptotically, if w G Wo, we have 



K-\w,b) 
These are reversed if w G W\. 



Re^-^ + OU-^) be B 



zlm-^-v + O 



7r(6— w) v 1 6— ui| 2 



6 G B x . 
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Figure 9: 



4.4 Bounded domains 

We can similarly compute the asymptotics of K~ x on bounded domains. 

Let U be a Jordan domain. Let L^e be a domain in 2eZ 2 which approximates U. It is 
a well-known fact that the discrete Dirichlet Green's function on L^e (defined by: y) 
is the expected number of passages at y of a simple random walk started at x and stopped 
at the boundary, alternatively, G(x,y) is the inverse of the discrete Laplacian with Dirichlet 
boundary conditions) approximates the continuous Dirichlet Green's function on U. 

We can use this fact to construct a graph U e G eZ 2 for which the K~ x operator converges 
nicely to the corresponding continuous object on U. 

The graph U e has a black vertex for each vertex and face of U^e, and a white vertex for 
each edge of U 2e (and U e is the superposition of C/ 2 e an d its dual). However we remove from 
U e one black vertex 6 on the outer face. 

Lemma 3 U e has the same number of black and white vertices. 

This follows simply from the Euler formula V — E+F = 1 applied to U^- In fact spanning 
trees of L^e rooted at bo are in bijection with perfect matchings of U e (see |T7J Ej). 

Let G(u, v) be the continuous Dirichlet Green's function on U. Let G(u, v) be the analytic 
function of v whose real part is G(u, v). On U e we have K^ 1 = 2edG + o(e), in the following 
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sense. Define F ,Fi and F + ,F_ by 



1 



dG(u,v) = -(F (u 1 v)du x + iFi(u,v)du 



V) 



2 

= -(F + (u,v)du + F^(u,v)du). 

Theorem 3 When \w — b\ is not 0(e) we have 

{eReF (w, b) + o(e) w e W and b e B 

eiImF (w 1 b) + o(e) w G W and b G B\ 

cReFi(w, b) + o(e) u> G W\ and b G B 

eiImFi(w, b) + o(e) w G Wi and 6 G Si. 

For example on M 2 we have t>) = — ^- log(v — u), so dG = ^- dUx +^u y anc j ^ _ ^ _ 
^^3^-. As another example, on the upper half plane we have G(u,v) = — ^log^5f, from 
which we get 

~ If du du 

d&(u, v) = — ■_ 

2n \v — u v — u j 

As a consequence F + (u, v) = n ^- u ) an d F_(u,v) = — ■ 

It is easier to work with the functions F± = F ± Fi since they transform nicely. Recall 
that G is conformally invariant, that is if <fi : V — > U is a conformal homeomorphism then 
Gy(u, v) = Gu(4>(u),4>(v)). As a consequence F± are conformally covariant: 

FX(u,v) = FK<P(u)A(v))<PXv). 



F v _{u,v) = F^((f)(u),(f)(v))(f)'(v). 

This allows us to compute asymptotics of K" 1 on any Jordan domain. 

More precisely, the result is the following. When w, b are far apart, we can use the above 
asymptotic formulas. When they are close together, the dominant term is controlled by (0) 
with a deviation of order e given by the function(s) F + (u,v) — n ^_ u ^ and F_(u, v). 

4.5 Height moment 

Here is a sample computation, the expected value of h(p)h(q) in the upper half plane. 

Let cij, 6j be the horizontal edges on a vertical path from p to the x-axis, the ctj being 
those with a white left vertex, and similarly let q, dj be those for q, see Figure [nil We have 
Hp) =^2~2i a i- &« h{q) = 4 J2j Cj - dj. So 

E(h{p)h(q)) = 16^E(a i c i ) - E(6 i c i ) - E(a i d i ) + Efadj). 

Each of these terms can be written as the determinant of a 2 x 2 matrix with entries in K~ x . 
For example 



E(ajCj) = det 
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Figure 10: 



where af \ (ip are the white and black vertices of and likewise for Cj. The diagonal terms 
will all cancel out, and after a certain amount of rearranging, one finds in the limit 



E((%) -%))(%)-%))) 



71,72 



F + ( Zl ,z 2 ) 
F + (z 2 , Zl ) 



dz\dz 2 - 



71:72 







F-( Zl ,z 2 ) 



7l>72 



F_(z u z 2 ) 
FJz 2 , Zl ) 



dz\dz 2 + 



F-(z 2 , Zl ) 
F + ( Zl ,z 2 ) 

'71,72 F +(Z2,Z!) 
- 1 T? („ „ \ _ 2 



dzidz 2 



dz\dz 2 . 



For the upper half-plane we have F+(zi, z 2 ) = 7r(za _ Zi) and z 2 ) = -^—=^. Plug- 

ging these in gives 



7T- 



1 ii 4 

dz\dz 2 + 



71 ^72 



(*2 " ^l) 2 



7T Z 



71 ^72 



(^2 - 



; d Zl dz 2 + 



4 

+— 

7T 2 



-71 ^ 72 1^2 - ^lj 

The first of these integrals gives 



dz\dz 2 — — ^ 



7T 



71 ^72 



(- 2 -zr/ ZldZ2 - 



4 log (p-^)( r - s ) 



7r 2 (p — s)(r — q) 



Therefore 



E((%) -%))(%)-%))) = -i f-2Relog {P S l + 2Relog {P ' °' )t/ " ^ 



7T Z 



(p- s)(r - g) 
P- V 



(p-s)(r - q) 



g 

—Re log 

Note that this quantity is a multiple of the Dirichlet Green's function G(p, q). 
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4.6 The Gaussian free field 



Let ho(p) = h(p) — h(p) be the fluctuation of h away from the mean. More generally one finds 
that all the moments of the height fluctuations can be written in terms of Green's functions: 

Theorem 4 Let U be a Jordan domain with smooth boundary. Letpi, . . . ,pk G U be distinct 
points. If k is odd we have lim e _ + oE(Ao(pi) • • • ^o(Pfc)) = 0. If k is even we have 



( 16\ fc/2 

UmE(/io(pi)---/io(p*)) = ( J G (Pt{i),P<t(2)) ■ ■ ■C(p (T(fc _i ) ,p (T(fc) ; 

pairings a 



By Wick's theorem, this implies that ho(x) is the unique Gaussian process with covari- 
ance function K(ho{x)ho(y)) = G{x,y). This can be taken as a definition of the massless 
Gaussian free field. 

An alternate description is as follows. 

On U let {ej} be an orthonormal eigenbasis of the Laplacian with Dirichlet boundary 
conditions, and Ae^ = A^. 
Define 



oo 

GFF(x) = J2^=ei{x) 



where c, are independent Gaussians of mean and variance 1. 

This series defines a distribution, not a function: the series diverges almost surely almost 
everywhere. However for a smooth function ip the series 



oo „ 

GFF(/0) d = ~7= / e^^dx 



converges almost surely. 

Theorem 5 ( [12J ) As e tends to 0, ho tends weakly in distribution to 4/a/7t times the 
massless free field GFF on U in the sense that for any smooth function ip on U, the random 
variable e 2 XLei/ e ~ip(x)h (x) tends in distribution to -^GFFijp). 

It is worth pointing out that, although values of h are integral (after adding a non- 
random smooth function h), h converges to a continuous object. 

There is a lot more that can be said about the distribution of the value of h and ho at 
a point, or the joint punctual distributions. It can be proved that the distribution of h at 

a point tends, when scaled by ^/log i, to a Gaussian with variance equal to some universal 
constant. 



5 FK percolation on critical planar graphs 

Up to now we have been working on subgraphs of Z 2 . For the dimer model on certain other 
regular graphs, such as the honeycomb graph, similar results can be obtained. However 
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for more general periodic planar graphs, the situation can be more complicated. In this 
section we discuss a different but related model, the FK-percolation model, on an interesting 
family of planar graphs the "isoradial" graphs. For this family of graphs we have a surprising 
property of the partition function and measures that they depend only on the local structure 
of the graph, not its long-range order. 

In the subsequent section we return to dimers but again on this same family of isoradial 
graphs. 

In this section we will show that the partition function of the random cluster model has 
a particular form for this special family of planar graphs (which includes the case Z 2 ). 

Recall the definition of the FK-percolation (random cluster) model. Let G be a graph. 
The space of configurations is X = {w : E — > {0, 1}} (each edge is open or closed). 

Let v : E — > (1, oo) be an assignment of weights to the edges, and q > a constant. 

We define a probability measure on X, 

m«o = \<t n ^ 

e open 

where c is the number of connected componets of open edges of w, the product is over open 
edges of w, and Z is a normalizing constant, called the partition function: 

Z =zZ ( l C II "(e)- 

w£X e open 

Note that when q = 1 we are reduced to the standard peroclation model. More generally 
when q is an integer the model is equivalent to the g-state Potts model PQ; in an appropriate 
limit q — > the model is equivalent to the spanning tree model. 

For further information on the model see [Zj. 



5.1 Duality 

If G is planar let G* be its planar dual. Let v* be weights on edges e* of G*, defined by 
v*(e*) = q/y{e). 

Lemma 4 For any configuration of edges of G with c components, k edges and r cycles, we 
have c — r + k = V , where V is the number of vertices of G. 

Then we have the following identity between Z and Z*, the dual partition function: 
Proposition 1 

Z* = Z ■ q E ~ V /W, 
where W is the product of all of the edge weights W = YleeE K e )- 

Proof: Suppose A is a configuration of open edges, with c components and weight 
q^lleeA ^(e). The dual configuration has edges E — A and d components. Its weight is 

/ n n dWi^rLw. 

e&E-A e£E-A V ; e&A 
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where W is the product of all edge weights of G. But by the Lemma, c + A — d = V, 
so the weight of a dual configuration is a constant q E ~ v /W times the weight of the primal 
configuration. □ 

In particular the duality is an isomorphism between the corresponding probability 
spaces. As a consequence we have p e + p e * = 1, where p e ,p e * are the edge probabilities. 

5.2 Y — A transformation 

We can sometimes transform G while preserving the measure /i. See Figure fTTl 




Figure 11: 

Given a triangle in G with vertices Vi,v 2 ,v 3 and edge weights e 23 = a, = b,ei 2 = c, 
and a 'V with the same vertices and edge weights A, B, C, (as in Figure [TTJ we have 



event 


weight in A 


weight in Y 


Vi,v 2 ,V3 connected 
v i , v 2 connected 
v i, t>3 connected 
v 2 , t>3 connected 
none connected 


(l + o)(l + 6)(l + c) 
(l + c)(l + a+b+ab ) 

(1 + 6)(1+ a+ " +ac ) 
(l + a)(l+ b+ f bc ) 

i 1 a+6+c 1 ab+ac+bc+abc 
q q 2 


q + A + B + C + AB + BC + AC + ABC 
q + A + B + C + AB + ac+bc+abc 

q + A + B + C + AC+ AB+BC+ABC 
q + A + B + C + BC+ AB + A ^+ABC 
q + A + B + C+ AB+AC+BC I ABC 



Here the events on the left refer to connections via open edges outside the triangle. The 
central column gives the additional (multiplicative) weight due to the possible open edges 
inside the triangle. The left column gives the additional (multiplicative) weight due to the 
possible open edges inside the Y. 

The Y — A transformation preserves the measures on condition that the weights in the 'A' 
column be proportional to the weights in the Y' column. This gives 4 polynomial equations 
for the weights a, b, c, A, B, C. 
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Lemma 5 These equations have a solution only if 



q + ah + ac + be + abc = 



(3) 



and 



q 2 -q(A + B + C) + ABC = 0. 



(4) 



In this case the solution is 



A 



a 



B 



Q 
b 



C 



c 



(5) 



For certain types of graphs/weights, these transformations take on a particularly nice 
form: 

5.3 Isoradial embeddings 

This notion is due to Duffin [S] and Mercat JH]. We say G has an isoradial embedding 

if G is drawn so that each face is a cyclic polygon, that is, is inscribable in a circle, and all 
the circles have radius 1. 



In this case we can embed G* isoradially as well, using the circle centers (at least as long 
as the circle centers are in the interior of the corresponding faces of G). Each edge has a 
rhombus around it whose vertices are the vertices of the edge and its dual fFigure fT^j) . 

Let 9 be the half-angle of the rhombus and define 




Figure 12: Isoradial embedding of a graph 




where r = cos 1 ( v /g/2). 
Then one may check that 
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Figure 13: 



Lemma 6 A Y — A transformation preserves isoradiality under these weights. 

This follows from plugging in the weights v{e) into Q, and using the fact that a triangle 
exists only when 9 a + d\, + 9 C = tt/2. 

5.4 Periodic graphs 

A planar graph G is periodic if translations in Z 2 are (weight-preserving) isomorphisms. In 
this case G n = G/nl? is a toroidal graph. Let Z n be the partition function on G n . we define 
Z = lim^oo Z]l^ Gn \ This is the partition function per site for the infinite periodic graph 
G. 

5.5 Main theorem 

Here is the main result of this section: 

Theorem 6 There exists a function F q : [0, f] — > R such that for any periodic, isoradial 
planar graph G, with edge weights v as above, the partition function per site is 



The proof doesn't tell us anything about the function F q . However F q was computed 
(non-rigorously) by Baxter [TJ (for the graph Z 2 , and hence for any isoradial graph). For 
q G (0,4) the answer is 




edges in a f. d. 
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Any information about this integral would be greatly appreciated. 
Proof sketch: Each Y <-> A changes the partition function by 

q 2 

Zy = Z^——. 

aoc 

Define a new-fangled "partition function" 

Zq^ 

It is easy to check that Zy = Z&. 

The idea is to use Y — A transformations to turn the graph into a graph with large blocks 
which are copies of big pieces of Z 2 . Because G is periodic, each bi- infinite rhombus chain 
is parallel to copies of itself. We'll rearrange the graph so as to move these chains adjacent 
to each other in clumps, so that the resulting rhombus tiling consists of large rhombi each 
tiled with N 2 copies of smaller versions of themselves (Figure 

How do we do this transformation? Convert rhombus chains into strings as in Figure ITU 



Figure 14: 

Slide the strings around without creating new intersections. When a string crosses the 
intersection of another pair of strings, the crossing corresponds to a Y — A transformation 
on the underlying graph. We can slide the strings around until the graph looks like a union 
of blocks, each tiled by parallel rhombi, as in Figure ITol That is, we can do this unless some 
rhombus chain is separated from its translates by a third chain which does not intersect 
either but whose rhombi have a different common parallel. In this case it is possible to show 
that simply exchanging the chains does not change the partition function. (One way to see 
this is to add an extra edge and it's "negative" at the same place so that you can exchange 
the two rows using a sequence of Y — A moves...) 
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Figure 15: 



Now by subadditivity, when there are only large blocks, the partition function can be 
computed separately on each block and the results multiplied together (with an error of 
lower order). But on each block the graph is a piece of 1? and the partition function is 
exp(Ac(l + o(l))) where A is the area and c a function of the angle only. Letting iV — > 00 
we can conclude in that Z = Hedges F q {9) f° r some Fq- The result follows. □ 

This theorem has some nontrivial consquences even without knowing F q . For example 
let Z rect be the partition function per site for the graph Z 2 with edge lengths 1 vertically and 

horizontally (and corresponding weights v as above). Let Zhex be the partition function 
per site for the honeycomb graph with edge lengths 1. Then 



J rect 



2cos(|cos- 1 (^)) 



As another "application" of the theory, it is natural to conjecture that the FK-percolation 
model is critical for weights v. For example in standard percolation (q = 1) we should have 

Pc/(1 -Pc) = V. 

For example on the on the honeycomb lattice the critical probability should be 

2 cos - 

p c = w 0.652704. 

yc 1 + 2 cos f 

This was in fact proved for the honeycomb lattice by Wierman j!9j . 

Proposition 2 The only function f for which (a, b, c) = (f(0 a ),f(9b),f(0 c )) satisfies (0) 
when 9 a + 0& + 9 C = tt/2 and (0) is 

sin (^9) 



where r = cos _1 ( v /g/2). 
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6 Integrability and dimers on critical planar graphs 



Surprisingly, the same isoradiality condition which worked so well in the random cluster 
model also has a simplifying effect on the dimer partition function. 

Let G be a bipartite planar graph. We embed G isoradially. Let v[e) = sin 29, where 9 
is the half-angle of the corresponding rhombus. We use these weights to define a probability 
measure on perfect matchings of (finite subgraphs of) G as we did before: the probability of 
a matching is proportional to the product of its edge weights. 

We define a Kasteleyn matrix as follows: if e = wb is an edge with dual edge e* = p*q*, 
then K(w, b) = i(p* — q*). 



Lemma 7 K is a Kasteleyn matrix for (any finite simply connected subpiece of) G. 

Note that |lif(tu,&)| = 2sin0. 
Proof: It suffices to show that K is Kasteleyn-flat [T^], that is, we must show that for 
each face of G\ with vertices u%, v i, . . . , u m , v m in cyclic order, we have 

arg(K(w 1; v{) . . . K(u m , v m )) = arg((-l) m_1 ir(w 1 , u 2 ) . . . K[v m -x, u m )K(v m , m)). 

(This identity implies that two dimer configurations which only differ around a single face 
have the same argument in the expansion of the determinant. By jT3], any two configurations 
can be obtained from one another by such displacements.) To prove this identity, note that 
it is true if the points are regularly spaced along a 2m-gon, and note that it remains true if 
you move one point at a time. □ 

There are two main results in this section, an explicit computation of Z and K . 

Theorem 7 The determinant per site of K satisfies 




Figure 16: 




edges e 
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where v[e) = 2sin(0(e)) ; N is the number of vertices in a fundamental domain, the sum is 
over the edges in a fundamental domain, and L is the Lobachevsky function, 

f ' 

L(x) — — \og2sintdt. 
Jo 

To describe K -1 , first define for any complex parameter z an elementary discrete 
analytic function to be a function satisfying f v * = f v /(z — e ta ), if vv' leads away from a 
white, or towards a black vertex, and f v i = f v ■ (z — e ia ), if vv' leads away from a black, 
or towards a white vertex. Such an / is well-defined the product of the multipliers going 
around a rhombus leads back to the starting value. 

Theorem 8 We have 

K^fawo) = -— ^ 6 ' ReS z=e^(fb/fw ), (6) 
poles e tS 

where the angles 9 G K are chosen appropriately. This can be written 

' ' log zdz, 



4n 2 i J c f W0 (z) 

where C is a closed contour surrounding cclw the part of the circle {e ld \ 9 e [9 Q — 7r + e,9o + 
7i — e]} which contains all the poles of f , and with the origin in its exterior. 

Proof: Let F(b) denote the right-hand side of ©• We will show that YlbeB K( w i b)F(b) = 
6 Wo (w), and F(b) tends to zero when b — > oo. 

Let Cg be a small loop around e td in C. We have 

k 

J2K(w,b)F(b) = £i(- e «J + e «J-i)F(6,), 
beB j=i 



1 - 

£i(_ e *j +c *j-i) ^ ^.Res z=e!9 (/ fe 



2tt? 

i=l poles e l£> 



i7r _ ~f 27TZ J Cfl 



2tt 

poles e l9 J=l 



2tt ^.„2vri 

poles e 



j _^ ^ /• / fc j 

~2^ 2m JcJ w ^ (z - e^-i) (:-«'"•) 

poles e«< 17 °9 \J=1 V 

L £ 



1 



9_ 

2tx ^ 2tvl 



poles e 



Odz = 0. 

C 9 
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However when w = u>o we have 

1 ( Oj- fy-i 



F{b 3 



2ni \ e i9 3 - e^- 1 
so that 

Yi{-&+&->.)F{b s ) = ±--£9,-0,-1 = 1, 

3=1 3=1 

since the angles increase by 2tt around wq. 

To complete the proof, one must show that F(b) — > as b — > oo. 

Proposition 3 We have 

K ~ 1 M = ^(t^— + 7T :L ^)+ 



2n \b — w ' (b — w) J ' \\b — w\ 2 J ' 
where 7 = /&(0) = e - *( 01+fl2 ^ rje^ - ^ J ' +a ^, one? i/ie angles are as illustrated in Figure [T9j 



u l 






u 4 






9/k// 










V u 3 / 










M a, 






a 4 

u 4 / 



w 



a. 



v 2 V, 



Figure 17: 



6.1 determinants 

For a finite graph G we define the normalized determinant, or determinant per sit 
deti M of an operator M: C G -> C G to be 

detiM^ldetMl 1 ^!, 

where IGI is the number of vertices of G. 
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For an operator M on a finite graph, detM is a function of the matrix entries M(i,j) 
which is linear in each entry separately. In particular for an edge e = ij, as a function 
of the matrix entry M(i,j) we have detM = a + (3M(i,j), where (3 is (— times the 
determinant of the minor obtained by removing row i and column j. That is, 

d(M(i,j)) 

or 

a(logdetM) 

<9M(z,j) Uj 

Suppose that G is periodic under translates by a lattice A. Let G n = G/nA, the finite 
graph which is the quotient of G by nA. Now if we sum (J7J) for all A-translates of edge ij in 
G n , and then divide both sides by \G n \, it yields 

dM{i,j) \d\ 

where M(i,j) is now the common weight of all translates of edge ij, that is, the left-hand 
side is the change in logdeti M when the weight of all translates of ij changes. Note that 
this equation is independent of n. 

It is now a short computation to compute the determinant of K from the exact form of 
IT -1 . ...see [13] 



6.2 Isoradial embeddings 

The set of isoradial embeddings ISO(G), when parametrised by the rhombus angles, is con- 
vex. This follows because an isoradial embedding is determined from a set of rhombus angles 
by the linear conditions that the sum around each vertex must be ir. 

There is a unique (possibly degenerate) isoradial embedding maximizing Z, because Z is 
strictly concave on ISO(G). 

(Joint with J-M. Schlenker) A zig-zag path in a planar graph is a path which turns 
maximally left at a vertex, then maximally right at the next vertex, then maximally left, 
and so on. That is, it leaves from a vertex along an edge which is adjacent in cyclic order 
to the edge it entered on, alteranting to the right and to the left. A planar graph has an 
isoradial embedding if and only if the following two conditions are satisfied (Figure EH) : 

1. No zig-zag path crosses itself 

2. No two zig-zag paths cross more than once. 

It is not hard to see why these conditions are necessary: a zig-zag path corresponds to a 
rhombus chain in an isoradial embedding, and such a chain is monotone and so cannot cross 
itself. Two such chains cannot cross more than once due to the orientation of their common 
rhombi. 



6.3 Hyperbolic ideal polyhedra 

To an isoradially embedded graph G we can associate an ideal polyhedron in H 3 as follows. 
Suppose G is embedded on the xy plane which is the boundary at oo of the upper half-space 
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Figure 18: Not allowed in isoradial graphs 



model of H 3 . Vertices of P are vertices of G* (circle centers) and edges of P are geodesies 
which project vertically to the edges of G*. 

The volume of P per fundamental domain can be seen to be the same as the entropy per 
fundamental domain of the dimer model on G. That is, 



can be associated with the mean curvature of P, which is by definition the sum over the 
edges of the dihedral angle times the (normalized) hyperbolic edge length. 

Question: Can one construct a tiling in some canonical (entropy-preserving) way using the 
geodesic flow inside PI 
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